Measurement of Propagation Time Dispersion 

in a Scintillator 



P. Achenbach*, C. Ayerbe Gayoso, J. Bernauer, R. Bohm, 
M.O. Distler, L. Doria, J. Friedrich, H. Merkel, U. Miiller, 
L. Nungesser, J. Pochodzalla, S. Sanchez Majos^, S. Schlimme, 
Th. Walcher, and M. Weinriefer 

Institut fiir Kernphysik, Johannes Gutenberg- Universitdt Mainz, Germany 



Abstract 

One contribution to the time resolution of a scintillation detector is the signal time 
spread due to path length variations of the detected photons from a point source. 
In an experimental study a rectangular scintillator was excited by means of a fast 
pulsed ultraviolet laser at different positions along its longitudinal axis. Timing 
measurements with a photomultiplier tube in a detection plane displaced from the 
scintillator end face showed a correlation between signal time and tube position 
indicating only a small distortion of photon angles during transmission. The data 
is in good agreement with a Monte Carlo simulation used to compute the average 
photon angle with respect to the detection plane and the average propagation time. 
Limitations on timing performance that arise from propagation time dispersion are 
expected for long and thin scintillators used in future particle identification systems. 
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1 Introduction 



There has been recent interest in the timing performance of long and thin 
plastic scintillation detectors with demands on the time resolution of 100 ps 
and below, e.g. for the future PANDA experiment at FAIR [1]. An existing 
experiment with a challanging scintillator system is COMPASS with the in- 
stallation of the 2.8 m long and 4 mm thin recoil detector under way [2]. This 
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paper shows that the detector performance required by these experiments is 
close to the fundamental limit that arises from the statistical fluctuations in 
the process of the signal generation. 

At PANDA , the identification of charged particles with momenta up to several 
GeV/c will be performed by the detection of internally reflected Cherenkov 
(DIRC) light [3|. Any measurement of the time-of-propagation (TOP) in the 
DIRC counter needs a precise time reference and can be combined with the 
time-of-flight (TOF) information from long scintillator slabs in front of the 
radiator barrel of the DIRC. Within the almost fully hermetic PANDA de- 
tector such a scintillator array needs to be as thin as physically possible. At 
COMPASS, it is the need for a low momentum proton detection that requires 
the scintillators to be very thin. 

The timing properties of scintillators are usually defined in terms of a coinci- 
dence time resolution. In a typical laboratory measurement the time difference 
between the discriminated signals of two photomultipliers (PMTs) placed at 
the extremes of a sample of scintillator is measured for minimum ionising par- 
ticles crossing the scintillator at its centre. For small counters of 45 cm length 
and 20 mm thickness a time resolution reduced for one PMT of the order of 
FWHM ^ 150 ps has been achieved by the authors in the laboratory. For long 
scintillators there are several effects degrading the time resolution, namely 
the light attenuation and the consequences of the scintillator acting as a light 
guide with corresponding propagation time dispersion. 

First photoelectron timing errors have been evaluated for slow scintillators 
(BGO, Tdecay ~ 300 iis) and larger amplitudes (number of photoelectrons 
N > 30) previously [I; but not for fast plastic scintillators and small 
amplitudes, where the achievable timing performance depends also on path 
length variations. The purpose of this paper is to discuss expressions for the 
timing error which include propagation time dispersion and to verify experi- 
mentally the the correlation between signal time and average photon angle at 
the read-out end of a scintillator. 



2 The Origin of Time Jitter in Scintillator Timing 



One fundamental limit on the time resolution of scintillation counters comes 
as a consequence of the statistical processes involved in the generation of the 
signal. Post and Schiff have first discussed such limitations. In general, 
the voltage pulse of a photomultiplier, Vrmt, can be expressed as a linear 
superposition of single photoelectron pulses, Vi, whose amplitude is allowed 
to vary due to gain fluctuations, gi, governed by a measurable probability 
distribution. The pulses arrive at individual times due to the time spread in 
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the energy transfer to the optical scintillator levels, tdep, the decay time of 
the light emitting states, temit, the propagation time, tprop, and the transit 
time, ^TT; being the time difference between photo-emission at the cathode 
and the arrival of the subsequent electric signal at the anode. In addition, 
there is white Gaussian electronic noise, w{t). We incorporate these processes 
into a general model: VpMT{t) = Eili QiViit - {Uep + temit + tprop + iTT)i) + 
w{t), where N fluctuates from one pulse to another. In a semi-classical model 
the probability for observing photoelectrons during a time interval T is 
given by a Poissonian distribution P{N, N) with N being the mean number 
of photoelectrons per pulse. 

The following calculation is focused on the limitations arising from the contri- 
bution of the path length variation. Central to the calculation is the relation 
between initial axial angle, 6, and propagation time: t = Ln/{ccos6), for a 
point like source of light placed at a distance L from the end face of a cylin- 
drical scintillator with refractive index in a medium of refractive index riext, 
where 6'max represents the complement of the total internal reflection angle, 
^max = arccosnext/''^- In the case of a rectangular shaped scintillator the re- 
lation is still valid as can be seen from the following consideration: in any 
reflection the change in the velocity vector of a particular photon takes place 
in the direction perpendicular to the face it hits. The component of this vec- 
tor along the longitudinal axis remains unchanged during the whole motion 
and therefore the propagation time only depends on the initial axial angle, 
distance to the PMT, and index of refraction. Consequently, a detector dis- 
placed by a distance d from the end face of the scintillator would link in one 
linear dimension the angle and propagation time with an accuracy given by 
ag ~ + ttx/d, where is the thickness of the scintillator and ax the 

aperture of the detector. 

For a quantitative description including extended light sources, light atten- 
uation and refraction a full photon tracking simulation with the correct ge- 
ometry of the set-up and a good model of scintillation light excitation and 
emission was needed. A general purpose Monte Carlo program f?\ simulating 
light propagation was used to compute the average photon angle with respect 
to the detection plane and the average propagation time. 

Experimentally, the quality of a scintillator surface limits the conservation of 
photon angles. Small geometrical inhomogeneities can have a strong impact 
on the transmission of photons and then the correlation between detection 
position and propagation time is heavily affected. There is no straight forward 
method for the estimation of the parameters necessary for an appropriate sim- 
ulation of the surface quality. In fact, a Monte Carlo refractive index matching 
technique was developed only very recently to determine these parameters 
It was therefore mandatory to perform an experiment with a sample of plastic 
scintillator in order to investigate empirically the mentioned correlation. 
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In an experiment the discriminators and any noise in the electronic circuits 
contribute to the time spread. Further, the response time of discriminators 
may shift with signal amplitude ("time walk"). Time walk effects are relevant 
if a large range of amplitudes is discriminated and can get corrected by various 
means, either in hardware or software. 

It is a standard in nuclear physics to use constant fraction discriminators 
(CFD) and leading edge discriminators (LED), the latter showing a charac- 
teristic time walk. In many cases constant fraction timing provides the best 
time resolution with scintillation counters. The situation is significantly dif- 
ferent when the number of photons is small. In this case the sampling of the 
arrival time distribution is poor and the signal shape is very variable. Fig. [1] 
shows the trace of a PMT voltage pulse for low light intensities as measured 
with a digital oscilloscope and the result of a Monte Carlo calculation for N = 
2 photoelectrons. The voltage pulse of a single photoelectron was modelled 
by Vspcit) = ^ct2e-*'/^7/o°°^^e-*'/^'dt, where g is the gain of the PMT, c is 
the conversion factor for charge to voltage and r has been chosen to be 4 ns 
in order to match the signal of the PMT used in the experimental set-up. A 
large distribution of pulse shapes reflecting the different arrival times of the 
individual photons is possible and the confounding effect of the overlapping 
single photon responses leads to problems with leading edge and peak detec- 
tors. Thus, the timing performance in the pile-up case degrades considerable 
from single photon timing. 



3 Calculation of Fundamental Limits from Propagation Time Dis- 
persion 

The spread of propagation times is given by Atprop = Ln/c{cos^^ ^^max — 1) for 
the minimum and maximum propagation times of meridional rays, tmin = Ln/c 
and tmax = Ln / {c cos O^s^x) . The probability density function of the prop- 
agation times of photons, dA^/(A^dt), produced by an event at t = 0, can 
be calculated from the angular distribution of photons inside the scintil- 
lator, dA^ = 27rdcos6'. It follows that dN/dt = -2% Ln/ {ct'^) and = 
//^;:^(dA^/dt)dt' = 27r(cos6' - 1). Finally, the probability density for the ar- 
rival times at the photocathode due to path length variations is Pprop(i) = 
dN/{Ndt) = Ln / {{cos e- l)ct^). 

The arrival time probability density function for a realistic scintillator can be 
folded in as follows: P{t) = I{t')P^,^^{t - t')Pi^{t - t')dt\ where PA(t) = 
exp —tc/ (nA) is survival probability depending on the bulk absorption coef- 
ficient, A, in the material. The distribution of decay times of light emitting 
states is included with a light pulse shape of the form: I{t) = [(e"*/"^^ — 
g-i/n ~) _ -p^^ _|_ i?e~*/'^^/r3]/(l -|- R), where Ti and T2 are two fast decay 
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constants, T3 corresponds to the decay time of the slow component, and R 
is the ratio of the slow to fast components. For the parameters of emission 
time distribution values of ri = 0.9 ns, T2 = 2.1ns, T3 = 14.2 ns, and R = 0.27 
were chosen. Comparable values are given in Ref. [9|. The probability density 
function of the transit time spread of electrons in the photomultiplier, Ptts is 
not relevant for the following discussion, since micro-channel plate PMTs pro- 
vide an alternative with excellent timing capabilities. Analytical descriptions 
of the probability density functions describing the time jitter of photoelectrons 
generated at the photocathode have been discussed before, e.g. in [lo| . 

The probability for a photoelectron to be detected in the interval between t^i^i 
and t is p{t) = // . dN/ {Ndt)dt' . Assuming a Poissonian distribution for the 
probability that a definite number of photoelectrons, iVthr, has been accumu- 
lated leads to the signal time distribution S{t) = iV^g^(^^^"^^) [p(t)]^""-i[l - 

which is to be normalised for getting the probability density func- 
tion. Calculated distributions for signal times with L = Im, n = 1.58, for an 
increasing number of photoelectrons in the pulse, = {5, 10, 25, 50} and 
Athr = 2, as well as for single photoelectron pulses (x 20) are shown in Fig. O 

The limit on time resolution arising from statistics of photon detection can now 
be calculated using the relations for the expectation values (tp) = Jl^.'^HSdt 
and (tp) = Jl^.^'H'^Sdt, and the variance of the signal time Var(tp) = {t^) — {tp)'^. 
The time resolution (FWHM) for L= 1 m, n= 1.58 as a function of A^ is shown 
in Fig. [3]for a range of thresholds, Athr= {2,5,10,25,50}. The curve for Athr= 
2 follows a FWHM cx 1/A^ dependence. 



4 Experimental Set-up 



The basic idea was to measure the average photon angle with respect to the 
detection plane at a distance d = 1 cm. A schematic representation of the 
arrangement of scintillator and detection plane is given in Fig. |H The trajec- 
tories of two photons leaving the scintillator in the same point are shown. A 
BC-408 scintillator from Bicron of dimensions 32 x 10 x 2000 mm^ was used. 
The scintillator is characterised by a decay time of Tdecay = 2.1ns, a small 
admixture of a slower decay time, and a refraction index oi n = 1.58. It was 
important for the conservation of photon angles that the faces of the scintilla- 
tor were not wrapped or covered. A minimum contact area of the bar at only 
two edges was accomplished by placing the scintillator along the longitudinal 
axis of a black metallic tube of the proper dimensions. The two tubes of 1 inch 
diameter were of type R4998 from Hamamatsu with a fast time response and a 
linear focused dynode structure of 10 stages, leading to a transit time spread 
of only (Ttts = 160 ps. At one end the reference PMT was fixed keeping con- 
tact with the scintillator by a suitable set of springs. At the other end a set 
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of two linear positioning stages permitted the scanning of the detection plane 
behind the scintillator, see Fig. [51 A thin plate with a 3 mm wide aperture was 
mounted in front of the PMT window to limit the acceptance. A ^°Sr source 
and a fast pulsed ultraviolet laser from Horiba Jobin Yvon were used for the 
measurements. The laser had a pulse duration FWHM < 200 ps. The signals 
were digitised by a charge integrating analogue-to-digital converter (LeCroy 
1885F, 50fC/count) and by a LeCroy time-to-digital converter (LeCroy 1875, 
25ps/count). The TDC digitisation corresponds to an accuracy of 25/vT2ps 
~7ps. 

A full scanning of the amplitudes was performed by means of filters placed in 
front of the movable PMT. The constant fraction discriminators showed little 
time walk {At/AQ < 0.7ps/pC). It was fitted and used for the correction of 
the TDC information. 



5 Discussion of Results 

A time resolution of FWHM ^ 130 — 260 ps was achieved for low amplitudes, 
see Fig. [6l With increasing amplitude the resolution improved until it levelled 
off at about FWHM ^ 80 ps, demonstrating that at high intensities the photon 
statistics is no longer decisive and the resolution is supposedly dominated by 
electronic noise. A function of the form FWHM oc + ai/{A + Aq) has been 
fitted to the data under the assumption that the detector response becomes 
Gaussian as the number of photoelectrons increases. 

Fig. [7] shows the positions of the aperture in the detection plane as points 
in the x- and y-dimension. For each point signal time and amplitude have 
been measured at a distance of 1 cm from the end face of the scintillator. 
Only the upper half plane was scanned due to the symmetric behaviour of all 
observables expected on grounds of the geometrical symmetry. The points are 
superimposed on the measured data presented as contour lines with steps of 
100 ps in signal time (right) and steps of —100 channels in amplitude (left). 
The signal time variation takes place earlier in the y-direction. The width of 
the scintillator in this direction is 10 mm to be compared to 32 mm for the 
perpendicular one. 

To exclude any amplitude dependent time shifts this experiment was repeated 
for three different positions of the laser source, adjusting the laser light in- 
tensity to equal PMT output signals for the central aperture position. Fig. [S] 
shows the variation of the signal time as a function of the position of the 
aperture with respect to the centre, when the laser light is injected at L = 
80, 100, and 140 cm. A signal time difference approximately linear with the 
distance for a given aperture position was observed. 
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In a perfect scintillator the portion of the light cone that lies inside the to- 
tal internal reflection angle is transported undistorted to the end of the slab. 
Since the transmission involves many reflections, the exit angle from a real 
scintillator is dependent on the sharpness of its edges and the parallelism of 
its faces. The measurements showed an unambiguous correlation between sig- 
nal time and tube position and demonstrated an sufficient conservation of 
axial angles. The observed correlation is not masked by a change of apparent 
anode sensitivity with angle of incidence at the photocathode. Although, the 
photomultiplier tube sensitivity depends on the incident angle since the path 
length of a photon through the photocathode increases like 1/ cos 6* giving 
a corresponding higher probability for the emission of a photoelectron 
The subsequent course of this electron is, however, not simple and the cosine 
law can only be considered as a first approximation. It would imply that at 
non-central aperture positions amplitude losses get partly compensated and 
effects of the propagation time dispersion get emphasised. The angular vari- 
ation of quantum efficiericy and the refiectance of a PMT entrance window 
were measured recently 1^, showing a larger photon-to-photoelectron con- 
version efficiency for incident angles 6 > 45° than that at ^ > 0°, however, the 
reflectance is also signiflcantly larger at these angles. 



Fig. [9] shows measured signal time (top) and amplitude (bottom) as a function 
of horizontal and vertical aperture position for the PMT at central position. 
Empty circles represent mirrored points. An analysis of the data reveals that 
the observed change in signal time is not following the change in signal ampli- 
tude. The light propagation is simulated in detail with the Monte Carlo code 
where first photoelectron timing was implemented. The observed variation of 
signal time and amplitude with position was verified, as shown in Fig. [TOl The 
code was then used to simulate the average photon angle with respect to the 
detection plane as a function of the horizontal and vertical aperture position, 
see Fig. [TTJ The good agreement between data and Monte Carlo together with 
the simulated relation between aperture position and angle verified that the 
propagation time dispersion has been observed. The measured time differences 
are in fair agreement with the spread of propagation times when applying the 
simple tA — ts = Ln/c{cos~^ 9a — cos~^ 9b) formula to a pair of simulated 
incident angles 9a/b at two aperture positions xa/b or i/a/b- 



The code was also used to simulate the time resolution of a full coverage PMT 
as a function of the number of photoelectrons. These simulations served as a 
guide for detector developments 13|. One example, where the contribution of 
the propagation time spread (Jprop gets significant is given by the cylindrical 
TOF counter being developed for PANDA [7]. A thickness of only 5 mm would 
allow to mount the scintillator strips together with the DIRC radiators, but 
reduces the number of photoelectrons to (A^pc) — 100 for particles crossing 
a scintillator at L = Im. The Monte Carlo simulation including the finite 
decay time of the scintillator predicts a minimum achievable time resolution 
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of FWHMjnin — 150 ps, depending a little on threshold. 

These results have brought forward the idea of improving the time resolution of 
a scintillation counter through position sensitive photon detection. An analogy 
of this method is used successfully in DIRC-like detectors using glass slabs for 
measuring the Cherenkov angle, but was never investigated in scintillators. 

Lower bounds on scintillation detector timing performance have been devel- 
oped in the past for moderate or relatively high total light output, addressing 



the exponential decay of the light intensity and non-ideal photodetectors [14 
These bounds can be useful in determining performance sensitivity to scin- 
tillator and PMT parameters. To the authors' knowledge the path length 
variation in long scintillators as a source of additional time dispersion has not 
been included in the models. 
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Fig. 1. Left: Typical PMT voltage pulse for low light intensities as measured with 
a digital oscilloscope. Right: Typical voltage pulse calculated by a Monte Carlo 
method for two single photoelectron pulses separated by 5 ns. 
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Fig. 2. Calculated distributions of signal times for L = 1 m, n = 1.58, for an increas- 
ing number of photoelectrons, N, at a fixed threshold of iVthr = 2: = 5 (solid 
line), N = 10 (dashed line), N = 25 (dot-dashed hne), and N = 50 (dotted line). 
The continuously falling (solid) line is a calculation for single photoelectron pulses 
(x 20). 
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Fig. 3. Calculated limits on time resolution (FWHM) for L = Im, n = 1.58, as 
a function of the number of photoelectrons, A^, at thresholds of A'thr = 2 (solid 
line following a FWHM cx 1/7V dependence), iVthr = 5 (dotted line), TVthr = 10 
(dot-dashed line), A'thr = 25 (dashed line), and A'thr = 50 (solid line). 




Fig. 4. Schematic representation of the arrangement of scintillator and detection 
plane at a distance d. The paths of two photons leaving the scintillator at the same 
point are shown. Typical parameters for the experiments were ax = a-y = 3 mm, 
L = 100 cm, d = 10 mm, = 32 mm, and ty = 10 mm. 
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Fig. 5. The photograph shows the crossed positioning units and the PMT equipped 
with a plastic mask with a square aperture of 3 mm x 3 mm. The PMT motion 
takes place in the plane perpendicular to the longitudinal axis of the scintillator. 




Fig. 6. Measured time resolution (FWHM) obtained with a BC-408 scintillator and 
an UV laser of varying primary light intensity as a function of the amplitude, A. 
The function FWHM oc uo + o-i/{A + Aq) has been fitted to the data. 
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Fig. 7. Measured signal time (right) and amplitude (left) as contour lines for the 
two dimensional scanning performed at a distance of 1 cm from the end face of the 
scintillator. The PMT positions for the data points are superimposed. The steps 
between contour lines correspond to 100 ps and —100 ADC channels, respectively. 
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Fig. 8. Measured signal time as a function of the aperture position when the scin- 
tillator is excited by the UV laser at three different positions, L = 80, 100, and 
140 cm. 
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Fig. 9. Measured signal time (top) and amplitude (bottom) as a function of hor- 
izontal and vertical aperture position. The width of the scintillator is 10 mm in 
y-direction and 32 mm in x-direction. Empty circles represent mirrored points. 
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Fig. 10. Simulated signal time (top) and amplitude (bottom) as a function of hori- 
zontal and vertical aperture position. The emission characteristics and the geometry 
of the scintillator have been adapted to the experiments. 
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Fig. 11. Simulated average photon angle with respect to the detection plane as a 
function of the horizontal and vertical aperture position. The emission characteris- 
tics and the geometry of the scintillator have been adapted to the experiments. 
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